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We study the vicinal surface of the restricted solid-on- 
solid model coupled with the Langmuir adsorbates which we 
regard as two-dimensional lattice gas without lateral interac- 
tion. The effect of the vapor pressure of the adsorbates in the 
environmental phase is taken into consideration through the 
chemical potential. We calculate the surface free energy /, 
the adsorption coverage O, the step tension 7, and the step 
stiffness 7 by the transfer matrix method combined with the 
density-matrix algorithm. Detailed step-density-dependence 
of / and Q is obtained. We draw the roughening transition 
curve in the plane of the temperature and the chemical po- 
tential of adsorbates. We find the multi-reentrant roughening 
transition accompanying the inverse roughening phenomena. 
We also find quasi-reentrant behavior in the step tension. 

PACS numbers: 05.70.Np, 68.35.Rh, 05.50.+q, 68.35.Md 



I. INTRODUCTION 

The adsorption effect has been an important sub- 
ject of study in wide area of surface science |[]-||. 
As for the steps on the surface, either with or with- 
out adsorption, detailed experiments in various length 
scales ftj|-|l2| has become possible due to the develop- 
ment of the surface observation technique such as STM 
(scanning-tunneling microscopy) |13|, LEEM (low-energy 
electron microscopy) |l4j and REM (reflection electron 
microscopy) [Gl|. The length scale of STM is the order 
of nanometer, while that of LEEM and REM is more 
than hundreds nanometer. Then, the importance of the 
detailed statistical mechanical study connecting observa- 
tions in different length scales has been increasing. 

As for the statistical mechanical study of vicinal sur- 
face with adsorption, however, reliable calculations have 
been difficult due to the criticality of the vicinal sur- 
face. In the present paper, to overcome the difficulty, 
we use the product wave-function renormalization group 
(PWFRG) algorithm p^Tq] which is a variant of the 
density matrix renormalization group (DMRG) algo- 
rithm p9| . The efficiency of the method for the vicinal 
surface problem has been demonstrated in Ref. ]16|,p7|,p0| 



We make statistical mechanical calculations on the vic- 
inal surface of the restricted solid-on-solid model coupled 
with the Langmuir adsorbates which we regard as two- 
dimensional lattice gas without lateral interaction. Wc 
aim at clarifying the "entropic effect" for surfaces and 
steps caused by adsorbates. The results are applicable to 
realistic vicinal surfaces if the temperature is high enough 
to neglect the interaction between the adsorbates, or the 
interaction itself is negligibly small. 

This paper is organized as follows. In §11, we introduce 
a "decorated" RSOS model where the ledge formation 
energy is locally modified by the presence of the adsor- 
bate atom on the ledge. After tracing out the lattice-gas 
degrees of freedom, we obtain the effective RSOS model 
which we show in §111. In §IV, we give phase diagram 
associated with the roughening transition. We show that 
the multi-reentrant roughening transition occurs. Calcu- 
lated results for surface free energy per projected area 
/, surface gradient p, magnetization (adsorption cover- 
age 9) M, step tension 7, and step stiffness 7 are given 
in §V. Quasi-reentrant behavior of step tension is also 
shown. In §VI, we give summary of the paper. 



II. MODEL HAMILTONIAN 

We consider the square lattice where the position of 
each site i is expressed in terms of a pair of integers as 
i = (m,n). To each site i we assign the integer height 
variable hi (— h(m,n)). An Ising spin variable a x (m,n) 
or a y (m,n) representing the adsorbate atom is assigned 
to each lattice edge (a x for vertical edge, and a y for hor- 
izontal edge). The Hamiltonian of the decorated RSOS 
model is written as follows: 

TC = ^^{e[l — aa y (m, n)] ■ \h(m + l,n) — him, n)\ 

+e[l — aa x (m, n)] ■ \h{m, n + 1) — h(m, n)\} 
— rfS^\h{m + 1, n) — h(m, n)} 



-Hy \<j x {m, n) + (T y (m, n)] 



(2.1) 
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TORSOS + + Wi, 



(2.2) 



where 



Wrsos = E - 

<»,i> 

— ?7 Vj[fe(m + l,n) — /i(m, n)], 

m,n 

= --ff [<j x (m, n) + a y (m, n)] , 

H m t = — J]] ea[a y (m, n) ■ \h{m + l,n) — /i(m, n)| 

+cr x (m, n) • |/i(m, n + 1) — /i(m, n)|]. (2.3) 

The Hamiltonian "Torsos describing the surface config- 
uration of the mother crystal is that of the restricted 
solid-on-solid (RSOS) model |2lJ] on the square lattice 
with nearest-neighbor (nn) interactions; e is the "bare" 
ledge energy, and r\ the Andr eev field [[22| to make sur- 
face tilt pG,B3|. In (2.z) and (2.3), the height difference 



hi — hj (= Ahij) is restricted to be A/iy = 0, ±1 (RSOS 
condition), which is a reasonable simplification because 
configurations with large |A/iy| are energetically unfa- 
vorable. 

The Hamiltonian Hl describes a system of free Ising 
spins (cr = ±1) on the edges of the square lattice (Fig. [I]). 
We regard these free spins as the Langmuir adsorbates in 
a monolayer adsorption without lateral interaction. The 
adsorption energy is represented by H. Since the adsor- 
bates are situated at the "bridge" sites, the total number 
of adsorption sites are twice as many as the total number 
of the lattice points of the RSOS model. 

The Hamiltonian 7i; n t represents the interaction be- 
tween the RSOS system and the adsorbate, with a being 
a parameter describing the local change of ledge energy 
caused by a spin (an adsorbate atom) (24). In what fol- 
lows, we assume a > — 1. 

The decorated RSOS model is a special case (J = 0) 
of the RSOS-Ising coupled model |2^,|2^] where the inter- 
action between adsorbates is taken into account. 

Since the lattice-gas representation is physically more 
natural for description of the adsorbates than the spin 
representation, we map the decorated RSOS to the equiv- 
alent RSOS-lattice-gas system. We introduce the lattice 
gas (LG) variable c = {0, 1} by c = (a + l)/2; a = 1 
corresponds to the presence of an adsorbate (c = 1), and 
a = — 1 to the absence of the adsorbate (c = 0). Then 
TL (eq. ( |2.2| )) is rewritten in terms of the LG variable as 
follows: 



Wlg = E e|fr-i ~ hj\ -rjy^\h(m + l,n) - h(m,n)] 

<i.j> m,n 

-/x ^ [c^m, n) + c y (m, n)] 

— ea[c y (m, n) ■ \h(m + l,n) — h(m, n)\ 



+c x {m, n) ■ \h(m, n + 1) — h(m, n)\] 
+const., 



where 



e = e(l + a), a 
\hi — hA = 



2a 
~ 1 + a 
{0,1}. 



2H, 



(2.4) 

(2.5) 
(2.6) 

The parameter e represents the microscopic ledge energy 
of the "clean" (i. e., without adsorbate) surface, fi the 
surface chemical po tent ial, and a the redefined coupling 
parameter. From ( Ej ) we see that, for a > (resp. 
a < 0), the adsorbate makes the microscopic ledge en- 
ergy smaller (resp. larger) than that of the clean surface, 
which will lead to rougher (resp. smoother) surface with 
larger (resp. smaller) population of adsorbates on the 
ledges (Fig. §. 

Physically, the chemical potential /i consists of two 
terms, 



A« = e a ds + k B T\nP, 



(2.7) 



where e a( j s is the adsorption energy which depends on 
the combination of the atomic species of the adsor- 
bate/substrate, and P is the environmental vapor pres- 
sure of adsorbate gas in thermal equilibrium with the 
surface {k B - Boltzmann constant, T: temperature). In- 
troducing Pq by e a ds = — k B ThiPo, we have 

H = k B Tln{P/P ) = 2H. (2.8) 

We should note the fol lowing symmetry of the sys- 
tem: the Hamiltonian ( |2.2|) is invariant under the change 
{a, cr x , c y , H) *-> {—ot, —<7 X , ~<^yi ~H). Although the lat- 
tice gas representation is more natural than the Ising-spin 
representation, the latter has a merit that this symmetry 
is apparently seen. In the following sections, we adopt 
the Ising-spin representation. Results in the lattice-gas 
representation are also shown if necessary. 



III. EFFECTIVE RSOS HAMILTONIAN 

Since each nn height difference Ah is restricted to 
{0, ±1}, we have an identity F(\Ah\) = F(0)5 Ahfi + 
F(l)(fiAh,i + <5a/i,-i) with Sij being the Kronecker's 
delta. Using this identity, we can exactly decimate spin 
variables {a} in H to obtain the effective Hamiltonian 
H eS as follows {fi = l/{k B T)): 



E 



exp{-0H) 



{er(m,n)} 



= (2cosh{3H) Nz / 2 exp(-(3H cS ), (3.1) 
(3H cS = Pe ee \h-hj\ 



<i,J> 



-77^[/i(m + l,n) -h{m,ri)], (3.2) 



,eff 



e — fee? 1 In 



cosh(/3ae + f3H) 
cosh (3H 



(3.3) 
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where N is the number of lattice point in the RSOS sys- 
tem and z (= 4) is the coordination number. Hence the 
decorated RSOS system is equivalent to a non-decorated 
RSOS system with the effective ledge energy e cS . 

In Fig. we show temperature and chemical poten- 
tial dependence of effective ledge energy e off . We should 



note that, at sufficiently low temperatures, 
approximated by 



,eff 



a - 



H/e\ + \H/e\). 



is well 



(3.4) 



We rewrite the e eS by using e and a in ( [2.5|) as follows: 



,eff 



e(l - a/2) - fc B Tln 



cosh(/3<5£/2 + (3H) 



cosh f3H 



At low temperatures, we have 



,eff 



e(l-a/2-\a/2 + H/e\ + \H/e\). 



(3.5) 



(3.6) 



In Fig. we show temperature and chemical potential 
dependence of e eff in the lattice-gas representation. 

At low temperatures where the entropy effect of adsor- 
bate configurations almost vanishes, e cff should bec ome 
approximately e(l — 59). Therefore, from eq. ( |3.6[ ), we 
have the low-temperature expression of O: 



9 = i + I(|a/2 
2 a 



H/~e\-\H/I\). 



(3.7) 



IV. 



MULTI-REENTRANT ROUGHENING 
TRANSITION 



(Fig. ^ (a)). In this limit, e eff becomes e(l — a) from 
which the roughening transition temperature becomes 
k B T R /e — (1 — oj)/Cr. In H — > — oo limit, we have 
(a) — »= —1, which corresponds to the clean surface 
(O = 0). In this limit, e off converges to e(l + a) — 
e, and the roughening transition temperature becomes 
k B T K /e = (1 + a)/Cn. 



B. Multi-reentrant transition 

Around a = 1, there occurs reentrant phase transition 
in a range of H (Fig. || (b) and (c)). At a = 1, as 
we decrease temperature fixing -ff/e = 0.01, the system 
cause phase transition twice: rough — > smooth — > rough 
(Fig. Q (b) and || (b)). This reentrant transition is similar 
to what has been known in the 2D frustrated Ising models 



At a = 1.1 and H/e = —0.0735, the phase transition 
occurs three times as we decrease temperature, which 
can be seen from Fig. ^ (c) and || (c); the system shows 
"multi-reentrant" roughening transition, which has not 
been known previously. The origin of this phenomenon 
is the non-trivial behavior of e eff in the parameter region 
around a = 1 and H = 0. In this region, e cS varies sensi- 
tively with H, making the critical curve non-monotonic. 

It should be noted that the (multi-) reentrant behav- 
ior accompanies the inverse roughening |2^] proposed for 
CsCl type crystal surface. Hence, our decorated RSOS 
model is a different but physically realizable system show- 
ing the inverse roughening. 



A. Determination of roughening transition 
temperature 

The roughening transition temperature of the original 
(non-decorated) RSOS model is known to be given by 
@ e/(k B Tg sos ) = 0.633... = Cr- Therefore, T R of 
the decorated RSOS model (the surface with Langmuir 
adsorption) is determined from e off /(^bIr) = Cr! when 
e cS /k B T > Cr (resp. < Cr), the surface is smooth (resp. 



(4.1) 



rough). From Eq. (|3.3|), we obtain an equation, 



__e f cosh[(a6 + H)/k B T R ] 

k B T R n \ cosh(iJ/fc B T R ) 



which determines Tr, (see Fig.0). For instance, at a = 0.5 



and H = 0, we obtain k B T R /e = 1.446. From Eq. ( |4.l| ) 
we have 



H 1 expjCR - (1 + «)e/(fc B T R )] - 1 
k B T R 2 n 1 - cx P [Cr - (1 - a)e/(k B T R )} ' 



(4.2) 



T 



which allows us to draw the critical line in the H 
plane, as shown in In Fig. |^. 

In H — > oo limit, we have (a) (magnetization) — >= 1; 
all of the ledge sites are covered by adsorbates (0 = 1) 



V. BEHAVIOR OF THE VICINAL SURFACE 
A. Universal behavior of the vicinal surfaces 

The exact mapping to the ordinary RSOS system given 
by ( |3.l| )^(3.3) implies that the vicinal surface of the 
decorated RSOS model belongs to the Gruber-Mullins- 
Pokrovsky-Talapov (GMPT) universality class |^-|3|] 
below the roughening temperature. This means that, in 
spite of the presence of the adsorbates, the coarse-grained 
vicinal surface is well described in terms of "effective ter- 
races" , "effective steps" and "effective kinks" allowing us 
to take the terrace-step-kink (TSK) picture. 

The surface free energy (per projected area) / = f{p) 
(p: step density) has the expansion characteristic to the 
GMPT universality class, 



f(p) = f(0)+7P + Bp 3 +O(p% 



(5.1) 



where 7 is the step tension and B is the step interaction 
coefficient. When the system is not isotropic, 7 and B 
depend on the mean running direction of steps. By 9, we 
denote the angle between one of the crystal axes on the 
facet plane and the mean running direction of steps. |36| ] 
We should then write 7 = 7(6*) and B = B(9). There 
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exists the following universal relation between 7(0) and 
B(9) |37],[58| (which leads to the universal Gaussian cur- 
vature jump at the facet edge for systems where 
the step-step interaction is short-ranged: 



B(0) 



tt 2 (fc B r) 2 

6 7(0) 



where 7(0) defined by 

7(0) = 1 {e)+d 2 1 (9)/d8 2 
is the step stiffness. 



(5.2) 



(5.3) 



B. 



Calculation of surface and step quantities by 
PWFRG method 



From the Hamiltonian (2.2) and the partition function 
Z associated with it, we obtain the Andreev surface free 
energy f(rf) as 



1 

'N 



^sTIii Z, 



(5.4) 



where if = (r] x , rj y ) is the Andreev field. As was shown by 
Andreev p2fl , f(ff) directly gives the equilibrium crystal 
shape (ECS). In the present study, we set r] x = rj and 
rj y = 0; y (77) — 77 curve can be regarded as the profile of 
ECS near (001) facet along the x-direction, and the facet 
size corresponds to the step tension 7(0) p9| . 

To calculate RHS of Eq.(5.4), we treat the system by 
the transfer-matrix method. Since the decorated RSOS 
model reduces to the original RSOS model, we can utilize 
the representation as the 19-state vertex model [^l],[l0| to 
make the transfer-matrix. We can also use the decorated 
RSOS model itself by mapping the model to a "deco- 
rated" vertex model (Fig. @) g||§ with 19 x 16 = 304 
non-zero vertex weights (304- vertex model). For approx- 
imate diagonalization of the transfer matrix, we employ 
the PWFRG method |l|-(l||2(J which is a variant of 



White's density matrix renormalization group (DMRG) 
method flf]. 

From f(ff), we can calculate f(p) (the surface free en- 
ergy with fixed surface gradient vector p), through the 
Legendre transformation: 



f(p) = f(v) + PV- 



(5.5) 



We should recall the following thermodynamical equa- 
tions associated with the Legendre transformation p2] |: 

Px = -df(ff)/dr] x , p y = -df(ff)/dr] y , (5.6) 
Vx = df(p)/dp x , % = df(p)/dp y . (5.7) 

In the actual calculation, we obtain p = p x {rj) as the 
statistical average 

p{rf) =< Ah x (m, n) >=< h(m + 1, n) — h{m, n) > . 

(5.8) 



By sweeping the field 77, we obtain ap-tj curv e. In orde r 
to obtain 7 and 7, we use the relations (5.1) and (fT^). 
Substituting Eq. (O) into Eq. (O), we have (for = 0) 



r] = 7(0) + 3£(0)|p| 2 + (higher order). 



(5.9) 



Hence, we obtain 7(0) and B(0) by performing the least- 
square fitting of the rj — p curve (inverted p — 77 curve). 
The step stiffness 7(0) is calculated from B(Q) through 
Eq. H). 



C. Surface free energies and surface gradient 

In Fig. H (b), we show f(rj) (r^ = 77, rj y — 0) at several 
temperatures for a = 0.5 and H = 0. The p — rj curves 
are shown in Fig. (7. The roughening temperature is 
determined by Eq. (O) as k B T R /e = 1.446 (§IV). At 



temperatures lower than Tr, the square-root behavior 
p ~ (77 — 7) 1 / 2 characteristic to the GMPT universality 
class is seen. 



From Eq. (5.5), we obtain f(p) which is shown in Fig. 
(a). At T < Tr, f(p) — p curve has a cusp at p = 0. 



D. Adsorbate concentration 



The magnetization M is calculated as M — (M x + 
My)/ 2 where M x =< a x > and M y —< a y > are di- 
rectly calculated by the PWFRG method. The adsorp- 
tion coverage O is given by O = (M + l)/2. We have the 
following expressions for M x and M y : 



M x = < <r x >= 



off! 



(5.10) 



ZMiMexPl-PH*] 
=<tanh/3i7f > h , 
My = < <j y >=< tanh f3Hy S (i, j) > h . 

Here 

Hf(i,j) = H + ea\Ah y (i,j)\, 

Ah y (i,j) = h(i,j + 1) - h(i,j), 

Hf(i,j) = H + ea\Ah x {i,j)\, 

Ah x (i,j)=h{i + l,j)-h(i,j), (5.11) 

can be interpreted as the effective local (magnetic) fields 
modulated by surface undulation. 

Since Ah takes or ±1, we have the identity < 
F(\Ah\) >= F(0) + (F(1)-F(0)) < \Ah\ >, which allows 
us to express M x and M y as 

M x = tanh/JH + [tanh/3(H + ea) 

-tanh/3/7] < \Ah y (i,j)\ > h) 
My = t&nh f3H + [tanh/3(i? + ea) 

-tanh/3i7] < \Ah x (i, j)\ > h ■ (5.12) 
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We should note the relations, 

< j)| >h = < £Afc B (i,j),i >fc 

< |A/i x (i, j)| > A = < <Ja/i„(»j),i >/i 

+ < ^Ah x (i,i),-1 >fci 



(5.13) 



where is the Kronecker's delta. Then, we have 

Py = < <W y (ij),l >h — < 8&h v (i,j)-X >h= 0, 

Px = < S A h x (i,j),i >h - < S A h x {i,j)-i >h= P- (5-14) 



Substituting Eq. and Eq. (p^) into Eq. ( p2| ) 

we obtain 

A4 = tanh/3# + [tanh /?(if + ecu) 
- tanh /3H] • 2 < 



M y = tanh/3# + [tanh /? (if + ea) 



tanh /?if] (p + 2 < (5 



A/t x (i,j), 



-1 >fc). 



(5.15) 



In Fig. we show p-dependence of M x and M y 
for a = 0.5 and H = 0. As is seen from the fig- 



ure, < 5 



Ah y {i,j) 



_l >h depends on p, taking maximum 



around p ~ 0.5. At p = 0, < <W„(i,j),-i ( =< 
>/j) is small because the surface is in the 



smooth phase; while at p 



1, < 6 



Ahy(i,j) 



-i >h= o 



due to the geometrical restriction in the RSOS model. 
Therefore, < 8Ah y (i,j),-i >h should have a maximum in 
between p = and p = 1 . 



E. Step tension and step stiffness 

At low temperatures, step quantities of the RSOS 
model are well approximated by the interface quantities 
of 2D Ising model pf|. Using the Ising- model results, we 
have the following approximate analytic expressions of 7 
and 7 El|,|l, 



rff 



(5.16) 



7 = e~ - fc B Tln[tanh(e ctt /2fc B T)], 
7 = fcBTsinh(7//cBT') 
= fcBT[exp(e off 7fc B T) tanh(e cff /(2fc B T)) 

- exp(-e cff /fc B T) coth(e cff 7(2fc B T))]/2. (5.17) 

In Fig. ^| andpp|, we show the temperature depen dence 
of st ep te nsion and step stiffness given by Eq. (5.16) 
and ( 5.17 ) for several values of H (or, equivalently, u/2). 
We see that the Ising- model expressions (|5.16|) and ( 5.17 ) 
agree very well with the PWFRG calculation at low tem- 
peratures. 

In the temperature dependence of 7, there is an inter- 
esting property that 7 exhibits the quasi-reentrant be- 
havior; there appears a temperature region where 7 in- 
creases as we raise the temperature, although 7 is usually 
a decreasing function of the temperature. Accordingly, 
7 has a maximum at a temperature above which 7 de- 
creases again as we raise the temperature. This behavior 



of 7 is what we call the quasi-reentrant behavior of a step, 
which can be seen in the 2D frustrated Ising models [E8| . 
Since the size of the facet in the ECS is proportional to 
the step tension, the facet size should also exhibit the 
quasi-reentrant behavior (28) . 

The step tension and step stiffness in the lattice-gas 
representation are given in Fig. [ll] and |l2|, respectively. 
When a = 0.5 (Fig. |l| (b) and Fig. y (b)), 7 and 
7 of the surface with adsorption are smaller than those 
of the clean surface; the step fluctuation width which is 
inversely proportional to 7 |4^j4j|] becomes larger, and 
Tr becomes lower. As the vapor pressure P is raised, 7 
and 7 become smaller because the increase of adsorbase 
concentration causes the decrease of e cff (Fig. || (b)). At 
a = —0.5, 7 and 7 of adsorbed surface are larger than 
those of clean surface, and Tr, becomes higher. As we 
raise P, 7 and 7 grows, because the increase of adsorbate 
concentration causes the increase of e cff (Fig. || (e)). 



VI. SUMMARY 

In this paper, we have discussed the Langmuir ad- 
sorption effect on the vicinal surface below the rough- 
ening temperature, in terms of the restricted solid-on- 
solid (RSOS) model coupled with the Ising spins. We 
have shown that this "decorated" RSOS model reduces 
to the ordinary RSOS model with the effective ledge en- 
ergy. The effective ledge energy depends on temperature 
T and the chemical potential fi of the adsorbates, which 
is caused by the change of the adsorbate concentration. 

We have found that, in a parameter region, there oc- 
curs novel multi-reentrant roughening transition. By em- 
ploying the product- wave- function renormalization group 
(PWFRG) method, we have calculated thermodynami- 
cal quantities associated with the vicinal surface of the 
model. We have also obtained quasi-reentrant behavior 
in the temperature dependence of step tension in a wide 
area of the phase diagram. 

We should give a comment on a related study made 
previously. In Ref. |44| , a similar modeling is adopted for 
the absolute SOS model and the discrete Gaussian model, 
where the adsorption effect on the roughening tempera- 
ture is discussed in terms of the approximate renormal- 
ization group approach (vicinal surface problem is not 
discussed). One of the merits of our RSOS modeling is 
that exact renormalization of coupling constant is avail- 
able, which allows us semi-analytical treatment of the 
problem. 
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FIG. 1. (a) RSOS heights [hi, . . . , hi) and the edge vari- 
ables in the mapped vertex model, (b) A vertex and leg vari- 
ables, (c) A vertex decorated by spins which are represented 
by circles. 

FIG. 2. Temperature dependence of the effective cou- 
pling constant e cS in the spin representation (Eq. (3.3)). (a) 
a — 0.1. (b) a = 0.5. (c) a = 1. In each figure, we draw 
the line for H/e = 0.2 (dashed), H/e = 0.1 (dotted), H/e = 
(solid), H/e = -0.1 (two-dot dashed), and H/e = -0.2 (dot 
dashed) . 



FIG. 3. Temperature dependence of the effective ledge en- 
ergy e eS in the lattice gas representation (Eq. (3.5)). In each 
figure, we draw the line for H/e — 0.2 (dashed), H/e — 0.1 
(dotted), H/e = (solid), H/e = -0.1 (two-dot dashed), and 
H/e = -0.2 (dot dashed), (a) a = 0.2. (b) a = 0.5. (c) 
5 = 1. (d) a = -0.2. (c) a = -0.5. (f) a = -1. 



FIG. 4. Determination of the roughening transition tem- 
perature T R (Eq. 4.1). (a) a = 0.5. H/e = 0.2 (dashed), 
H/e = 0.1 (dotted), H/e = (thick solid), H/e = -0.1 
(two-dot dashed), and H/e = —0.2 (dot dashed). (b) 
a = 1.0. H/e = 0.02 (dashed), H/e = 0.01 (dotted), 
H/e = (thick solid), H/e = -0.01 (two-dot dashed), and 
H/e = -0.02 (dot dashed), (c) a = 1.1 H/e = -0.0685 
(dashed), H/e = -0.0735 (dotted), and H/e = -0.0785 
(dot dashed). The horizontal line (thin solid) corresponds 
to Cr = e cB /k B T = 0.633. 



FIG. 5. Phase diagram, (a) a = 0.5. (b) a = 1.0. (c) 
a = 1.1. (d) a = 1.2. 



FIG. 6. (a) Temperature dependence of the f[p)/ksT. 
(b) Temperature dependence of the Andreev free energy f[rj) 
for a — 0.5 and H = 0. From top to bottom, ksT/e — 0.4, 
0.6, 0.8, 1.0 and 1.2. 



FIG. 7. PWFRG results for the p-r) curves (a = 0.5, and 
H = 0). Temperature of each curve is k^T/e = 0.4, 0.6, 0.8, 
1, 1.2, and 1.4, from right to left. 



FIG. 8. p-dependence of M x (solid) and M y (dashed) 
for a — 0.5 and H = 0. From top to bottom at p = 0, 
k B T/e = 1.1, 0.8 and 0.5. 



FIG. 9. Temperature-dependence of the step tension 7 
(Eq. (5.16)). Open circles correspond to PWFRG calculation 
for a — 0.5 and H = 0. We draw the line for H/e = 0.2 
(dashed), H/e = 0.1 (dotted), H/e = (solid), H/e = -0.1 
(two-dot dashed), and H/e = —0.2 (dot dashed), (a) a = 0.1. 
(b) a = 0.5. 



FIG. 10. Temperature-dependence of the step stiffness 7 
(Eq. (5.17)). Open circles correspond to PWFRG calculation 
for a = 0.5 and H = 0. We draw the line for H/e = 0.2 
(dashed), H/e = 0.1 (dotted), H/e = (solid), H/e = -0.1 
(two-dot dashed), and H/e = —0.2 (dot dashed), (a) a = 0.1. 
(b) a = 0.5. 



FIG. 11. Temperature-dependence of the step tension 
7 (Eq. (5.16) and (3.5)) in the lattice gas representation. 
We draw the line for H/e = 0.2 (dashed), H/e = 0.1 (dot- 
ted), H/e = (solid), H/e = -0.1 (two-dot dashed), and 
H/e = —0.2 (dot dashed), (a) The upper group of the line 
corresponds to a = —0.2 and the lower group to a = 0.2. (b) 
The upper group of the line corresponds to a = —0.5 and the 
lower group to a = 0.5. 



FIG. 12. Temperature-dependence of the step stiffness 
7 (Eq. (5.17) and (3.5)) in the lattice gas representation. 
We draw the line for H/e = 0.2 (dashed), H/e = 0.1 (dot- 
ted), H/e = (solid), H/e = -0.1 (two-dot dashed), and 
H/e = —0.2 (dot dashed), (a) The upper group of the line 
corresponds to a = —0.2 and the lower group to a = 0.2. (b) 
The upper group of the line corresponds to a = —0.5 and the 
lower group to a = 0.5. 
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